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(N '. Abstract 

l,_i'' Plutonium dioxide is of high technological importance in nuclear fuel cycle and is particularly 

o : 

"^ ' crucial in long-term storage of Pu-based radioactive waste. Using first-principles density-functional 

• , properties, and pressure induced structural transition of Pu02- To properly describe the strong 

Ch ' correlation in Pu 5/ electrons, the local density approximation+[/ and the generalized gradient 

'^ ■ approximation+[/ theoretical formalisms have been employed. We optimize U parameter in cal- 

O . culating the total energy, lattice parameters, and bulk modulus at nonmagnetic, ferromagnetic, 

and antiferromagnetic configurations for both ground state fluorite structure and high pressure 

(N ■ 

^ ■ cotunnite structure. Best agreement with experiments is obtained by tuning the effective Hubbard 

C^^ . parameter U at around 4 eV within LDA+f/ approach. After carefully testing the validity of the 

(N ; 

^^ \ ground-state calculation, we further investigate the bonding nature, elastic constants, various mod- 

vn ■ 

O ' uli, Debye temperature, hardness, ideal tensile strength, and phonon dispersion for fluorite Pu02. 

o: 

^— I ■ Some thermodynamic properties, e.g., Gibbs free energy, volume thermal expansion, and specific 

> ■ 

• ^H , heat, are also calculated. As for cotunnite phase, besides elastic constants, various moduli, and 

\-^ ' Debye temperature at GPa, we have further presented our calculated electronic, structural, and 



C^ 



magnetic properties for Pu02 under pressure up to 280 GPa. A metallic transition at around 133 
GPa and an isostructural transition in pressure range of 75-133 GPa are predicted. Additionally, 
as an illustration on the valency trend and subsequent effect on the mechanical properties, the 
calculated results for other actinide metal dioxides (Th02, UO2, and Np02) are also presented. 

PACS numbers: 71.27.+a, 61.50.Ks, 62.20.-x, 63.20.dk 



I. INTRODUCTION 



Actinide elements and compounds possess particular interesting physical behaviors due 
to the 5/ states and have always attracted extensive attention because of their importance 
in nuclear fuel cycle \Mi^- Actinide dioxides, AO2 {A=Th, U, Np or Pu), are universally 
used as advanced fuel materials for nuclear reactors and PUO2 also plays an important role 
in the plutonium re-use, separation and/or long-term storage. Recently, there has occurred 
in the literature a series of experimental reports |5|-[7| on the strategies of storage of Pu- 
based waste. Exposure to air and moisture, metallic plutonium surface easily oxidizes to 
PU2O3 and PUO2. Although the existence of Pu02+x {x <0.27) was reported by Haschke 
et al. [5|, recent photoemission study found that PUO2 was only covered by a chemisorbed 
layer of oxygen and can be easily desorbed at elevated temperature [7|, therefore, PUO2 is 
the stablest plutonium oxide. At ambient condition PUO2 crystallizes in a fiuorite structure 
with space group FmSm (No. 225). Its cubic unit cell is composed of four PUO2 formula 
units with plutonium atoms and oxygen atoms in 4a and 8c sites, respectively. By using 
the energy dispersive x-ray diffraction method, Dancausse et al |8| reported that at 39 GPa, 
PUO2 undergoes a phase transition to an orthorhombic structure of cotunnite type with 
space group Pnma (No. 62). 

As for the electronic-structure study of PUO2, conventional density functional theory 
(DFT) that applies the local density approximation (LDA) or generalized gradient approx- 
imation (GGA) underestimates the strong on-site Coulomb repulsion of the 5/ electrons 
and, consequently, describes PUO2 as incorrect ferromagnetic (FM) conductor |9| instead of 
antiferromagnetic (AFM) Mott insulator reported by experiment 10|]. Same problems have 
been confirmed in studying other correlated materials within the pure LDA/GGA schemes. 



Fortunately, several approaches, the LDA/GGA+f/ |lll4l3l|. the hybrid density functional 
of HSE 14|, and the self-interaction corrected local spin-density (SIC-LSD) |15|, have been 
developed to correct these failures in calculations of actinide compounds. The effective mod- 
ification of pure DFT by LDA/GGA+[/ formalisms has been confirmed widely in study of 



UO, 



11 



13| and PuOa 



16 



201 1 . By tuning the effective Hubbard parameter in a reasonable 



range, the AFM Mott insulator feature was correctly calculated and the structural parame- 
ters as well as the electronic structure are well in accord with experiments. Lattice dynamical 
properties of UO2 and PUO2 and various contributions to their thermal conductivities were 



studied by using a combination of LDA and Dynamical Mean-Field Theory (DMFT) |21 |. 
However, despite that abundant researches on the structural, electronic, optical, thermo- 
dynamic properties of plutonium dioxide have been performed, relatively little is known 
regarding its chemical bonding, mechanical properties, and phonon dispersion. In addition, 
the pressure induced structural transition has also not yet been studied by first-principles 
DFT-I- U calculations. 

In this work, we have systematically calculated the ground-state structural parameters, 
electronic, mechanical, thermodynamic properties, and pressure induced structural tran- 
sition of Pu02 by employing the LDA+ U and GGA+ U schemes due to Dudarev et al. 



11 



13| . The validity of the ground-state calculation is carefully tested. Our calculated 

attice parameter and bulk modulus B are well consistent with previous LDA-I- U results 

181]. The total energy, structural parameters, bulk modulus 5, and pressure derivative of 

the bulk modulus B' for AFM and FM phases calculated in wide range of effective Hubbard 

U parameter are presented and our electronic spectrum reproduce briefly the main feature 



of our previous study IGj. In addition, the bonding nature of A—0 bond in Pu02, Np02, 



UO2, and Th02 involving its mixed ionic/covalent character is investigated by performing 



the Bader analysis 22|, |23|. We find that about 2.40, 2.48, 2.56, and 2.66 electrons trans- 
fer from each Pu, Np, U or Th atom to O atom, respectively. In study of the mechanical 
properties, we first calculate the elastic constants of both Fm3m and Pnma phases and then 
the elastic moduli, Poisson's ratio, and Debye temperature are deduced from the calculated 
elastic constants. Hardness and ideal tensile strength of FmSm PUO2 are also obtained by 
LDA-I- U approach with one typical value of effective U parameter. The hardness of PUO2 
is equal to 26.6 GPa and the ideal tensile strengthes are calculated to be 81.2, 28.3, and 
16.8 GPa for pulling in the [001], [110], and [111] directions, respectively. As for the ther- 
modynamic study, the phonon dispersion illustrates the stability of PUO2 and we further 
predict the lattice vibration energy, thermal expansion, and specific heat by utilizing the 
quasiharmonic approximation based on the first-principles phonon density of state (DOS). 
One more aim of the present work is to extend the description of PUO2 to high pressures. 
The electronic, structural, and magnetic behavior of PUO2 under pressure up to 280 GPa 
have been evaluated. Results show that there occurs a metallic transition at around 133 
GPa for Pnma phase. The isostructural transition, similar to UO2 [24] and Th02 |25|], in 
pressure domain of 75-133 GPa is predicted. The rest of this paper is arranged as follows. 



In Sec. II the computational method is briefly described. In Sec. Ill we present and discuss 
our results. In Sec. IV we summarize the conclusions of this work. 

II. COMPUTATIONAL METHODS 



The DFT calculations are performed on the basis of the frozen-core projected augmented 



wave (PAW) method of Blochl 



27| using the LDA and GGA 



26 
281 



encoded in Vienna ab initio simulation package (VASP) 



29|. For the plane-wave set, a cutoff energy of 500 eV 



is used. The A;-point meshes in the full wedge of the Brillouin zone (BZ) are sampled by 
9x9x9 and 9x15x9 grids according to the Monkhorst-Pack 30|] scheme for fluorite and 
cotunnite Pu02, respectively, and all atoms are fully relaxed until the Hellmann-Feynman 
forces become less than 0.02 eV/A. The plutonium 6s'^7s'^6p^6d'^5f'^ and the oxygen 2s^2j»^ 
electrons are treated as valence electrons. The strong on-site Coulomb repulsion among the 
localized Pu 5/ electrons is described by using the LDA/GGA+ U formalisms formulated by 



Dudarev et at 



.0 



13j , where the double counting correction has already been included. As 



concluded in some previous studies of actinide dioxides (^02), although the pure LDA and 
GGA fail to depict the electronic structure, especially the insulating nature and the occupied- 



state character of UO-; 



M l3l|, NpOa \m, la^, and PuGs 



16 



17|, ll9|, the LDA/GGA+[/ 



approaches will capture the Mott insulating properties of the strongly correlated U 5/, Np 
5/, and Pu 5/ electrons in AO2 adequately. In this paper the Goulomb U is treated as a 
variable, while the exchange energy is set to be a constant J=0.75 eV. This value of J is 
same with our previous study of plutonium oxides [ly, [l7(]. Since only the difference between 
U and J is signiflcant [1^, thus we will henceforth label them as one single parameter, for 
simplicity labeled as U, while keeping in mind that the non-zero J has been used during 
calculations. 

Both spin-unpolarized and spin-polarized calculations are performed in this study. Com- 
pared to FM and AFM phases, the nonmagnetic (NM) phase is not energetically favorable 
both in the LDA+ U and GGA+ U formalisms. Therefore, the results of NM are not pre- 
sented in the following. The dependence of the total energy (per formula unit at respective 
optimum geometries) on U for both FM and AFM phases within the LDA+ U and GGA+ U 
formalisms is shown in Fig. [H At U=0 and 1.0 eV, the total energy of the FM phase is lower 
than that of the AFM phase either in LDA+ U scheme or GGA+ U scheme. However, as 
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FIG. 1: (Color online) Dependence of the total energies (per formula unit) on U for FM and AFM 

Pu02. 

shown in Fig. 1, it is clear that the total energy of the AFM phase decreases to become lower 
than that of the FM phase when increasing U. The total-energy differences {Epu—E afm) 
within the LDA+ U and GGA+ U at f/=4 eV are 0.705 and 0.651 eV, respectively. Both FM 
and AFM results will be presented in the following analysis. Besides, while the spin-orbit 
coupling (SOC) is important for certain properties of heavy metal compounds, it has been 
numerically found |9|, |33| and physically analyzed IJ, |3J] that the inclusion of the SOC has 
little effect on the bulk and one-electron properties of UO2 and PUO2. Our test calculations 
also show that within LDA+ U approach with 11=4 eV, inclusion of SOC will increase the 
optimum lattice constant by only 0.7% and the bulk modulus by about 0.5 GPa. Therefore, 
in our following calculations of plutonium dioxide, the SOC is not included. 

In present work, the theoretical equilibrium volume, bulk modulus B, and pressure deriva- 
tive of the bulk modulus B' are obtained by fitting the energy- volume data with the third- 
order Birch-Murnaghan equation of state (EOS) [351]. In order to calculate elastic constants, 
a small strain is applied onto the structure. For small strain e, Hooke's law is valid and the 
crystal energy E(y, e) can be expanded as a Taylor series |36|], 



E{V, e) = E{Vo, 0) + K) X^ ^.e, + y X^ Q,e,e, + 0{{e^}), 



(1) 



where E{Vq, 0) is the energy of the unstrained system with the equilibrium volume Vq, e is 



strain tensor which has matrix elements Eij {i,j=l, 2, and 3) defined by 

/ p ip ip \ 



2^6 ^2 2^4 
\ 2^5 2^4 63 



(2) 



/ 



and Cij are the elastic constants. For cubic structures, there are three independent elastic 
constants (Cii, C12, and C44). So, the elastic constants for fiuorite PUO2 can be calculated 
from three different strains listed as follows: 



{6, 6, 6, 0, 0, 0), e^ = {6, 0, 6, 0, 0, 0), e^ = (0, 0, 0, 6, 6, 6). 



(3) 



The strain amplitude S is varied in steps of 0.006 from 5=— 0.036 to 0.036 and the total 
energies E(y, 6) at these strain steps are calculated. After obtaining elastic constants, we 
can calculate bulk and shear moduli from the Voigt-Reuss-Hill (VRH) approximations 37 |. 
The Voigt (Reuss) bounds on the bulk modulus By (Br) and shear modulus Gy ( Gr) for 
this cubic crystal system are deduced from the formulae of elastic moduli in Ref. 38|. As 



for cotunnite PUO2, the nine independent elastic constants (Cn, C12, C13, C22, C'23, C'ss, 
C44, C55, and Cge) can be obtained from nine different strains listed in the following: 

e^ = ((5, 0, 0, 0, 0, 0), e^ = (0, S, 0, 0, 0, 0), e^ = (0, 0, S, 0, 0, 0), 



e^ = (0, 0, 0, S, 0, 0), e^ = (0, 0, 0, 0, S, 0),e^ = (0, 0, 0, 0, 0, 5), 
e^ = {5, 5, 0, 0, 0, 0), e^ = (0, 5, 5, 0, 0, 0), e^ = (5, 0, 5, 0, 0, 0) 



(4) 



are from Ref. [39|. Based 



and the formulae of elastic moduli in VRH approximations 

on the Hill approximation [37|, B=^{Br + By) and G=^{Gji + Gy). The Young's modulus 

E and Poisson's ratio v are given by the following formulae: 



E = 9BG/{3B + G),v = i3B- 2G)/[2{3B + G)]. 



(5) 



In addition, the elastic properties of a solid can also be related to thermodynamical pa- 
rameters especially specific heat, thermal expansion, Debye temperature, melting point, and 
Griineisen parameter 40|. From this point of view, Debye temperature is one of fundamen- 
tal parameters for solid materials. Due to the fact that at low temperatures the vibrational 
excitations arise solely from acoustic vibrations, therefore, the Debye temperature calculated 



from elastic constants is the same as that determined from specific heat measurements. The 
relation between the Debye temperature {6d) and the average sound wave velocity {vm) is 

where h and ks are Planck and Boltzmann constants, respectively, n is the number of atoms 
in the molecule and Vt is molecular volume. The average wave velocity in the polycrystalline 
materials is approximately given as 

-1/3 



Vr, 



1/2 



3 V vl vf 



(7) 



where Vt=\/G j p (p is the density) and vi=y^{3B + 4G)/3p are the transverse and longitu- 
dinal elastic wave velocities of the polycrystalline materials, respectively. 

III. RESULTS 

A. Atomic and electronic structures of fluorite Pu02 

We report in TableUthe lattice parameter (ao), bulk modulus (B), and pressure derivative 
of the bulk modulus {B') for AFM and FM Pu02 obtained in LDA+ U and GGA+ U frame- 



works. All these va^ 



B 



ues are determined by EOS fitting. For comparison, the experimental 



values of Oq (Ref. [5|]) and B (Ref. [4l|) are also listed. In the overall view, the dependence 
of the lattice parameter oq of AFM Pu02 on U is well consistent with our previous study 
[16| either in LDA+ U scheme or in GGA+ U scheme. For the LDA/GGA+ U approaches, 
the calculated oq improves upon the pure LDA/GGA by steadily increasing its amplitude 
with U. Actually, at a typical value U=4 eV, the LDA+U gives ao=5.362 A, which is very 
close to the experimental value of 5.398 A, and the GGA+U gives ao=5.466 A. Note that 

n n 

recent PBE+U [lO^ and LDA+U [18| calculations with U=4.0 eV and J=0.7 eV predicted 
the lattice parameter of AFM Pu02 to be 5.444 and 5.354 A, respectively, and the HSE 

n n » 

|14l | and SIC-LSD [15| calculations gave the values to be 5.385 and 5.44 A, respectively. On 
the other hand, as shown in Table [H the tendency of ao with U for FM Pu02 is similar to 
that of the AFM phase. At a typical value f/=4 eV, the LDA/GGA+[/ give ao=5.338 and 
5.452 A. These values are slightly smaller than those of the AFM phase. Previous HSE [ij] 
calculation of FM Pu02 gave the lattice parameter to be 5.387 A. As for the dependence of 



TABLE I: Calculated lattice parameters {ao), bulk modulus (B), and pressure derivative of the 

bulk modulus {B ) for AFM and FM Pu02 at GPa. For comparison, experimental values are 
also listed. 

Magnetism Method property U=0 U=l U=2 U=3 U=4: U=5 U=6 Expt. 

AFM LDA+f/ ao (A) 5.275 5.313 5.338 5.351 5.362 5.371 5.378 5.398'^ 

B (GPa) 218 208 224 224 225 226 227 178^* 

B' 4.1 3.7 4.3 4.3 4.3 4.3 4.3 

GGA+U ao (A) 5.396 5.433 5.446 5.457 5.466 5.473 5.480 

B (GPa) 185 188 191 192 193 194 195 

B' 4.3 3.8 4.4 4.5 4.5 4.5 4.4 

FM LDA+U ao (A) 5.270 5.290 5.309 5.325 5.338 5.350 5.361 

B (GPa) 230 224 221 220 218 215 212 

B' 4.4 4.4 4.3 4.4 4.4 4.4 4.5 

GGA+f7 ao (A) 5.384 5.405 5.424 5.439 5.452 5.464 5.476 

B (GPa) 193 188 186 184 182 179 174 

B' 4.5 4.4 4.5 4.5 4.7 4.7 4.8 



"' Reference 



Reference j4ll] 



bulk modulus B of AFM and FM Pu02 on U, it is clear that the LDA results are always 
higher than the GGA results, which is due to the overbinding effect of the LDA approach. 
While the bulk modulus B of AFM phase increases steadily with increasing the amplitude of 
f/, the bulk modulus B of FM phase will decrease along with the increasing of Hubbard U. 
At a typical value U=A eV, the LDA+ U and GGA+ U give 5=225 and 193 GPa for AFM 
phase, respectively, and 5=218 and 182 GPa for FM phase, respectively. Apparently, the 
GGA+[/ approach gives more close values to the experimental data of 5=178 GPa [41]. In 
our present study, as shown in Table [TTl the bulk modulus B deduced from elastic constants 



turns out to be very close to that obtained by EOS fitting. This indicates that our ca^ 
lations are consistent and reliable. In addition, as listed in Table Ull recent PBE+ U 



cu- 

and LDA+[/ llSll calculations with ^=4.0 eV and J=0.7 eV predicted the bulk modulus 
B of AFM PuOs to be 199 and 226 GPa, respectively, and the HSE [14] and SIC-LSD [l5| 
calculations gave the the bulk modulus to be 221 and 214 GPa, respectively. As for pressure 



TABLE II: Calculated elastic constants, various moduli, Poisson's ratio (v), spin moments {fJ-magX 
and energy band gap (Eg) for FmSm AFM Pu02 at GPa. For comparison, experimental values 
and other theoretical results are also listed. 



Method Cii Cu Cm B G E 

(GPa) (GPa) (GPa) (GPa) (GPa) (GPa) 



f^mag. ^g 

(mb) (eV) 



LDA+ U ( U= 


=0) 


386.6 


136.5 


71.9 


220 


89.9 


237.3 


0.320 


3.937 


0.0 


LDA+ U ( U= 


=4) 


319.6 


177.8 


74.5 


225 


73.0 


197.7 


0.354 


4.126 


1.5 


GGA+U {U= 


=0) 


343.9 


112.3 


53.7 


190 


73.5 


195.1 


0.328 


4.085 


0.0 


GGA+ U ( U= 


=4) 


256.5 


167.9 


59.2 


197 


52.7 


145.2 


0.377 


4.165 


1.5 


Expt. 










178^^ 










1.8^ 


FBE+U" 










199 








3.89 


2.2 


LDA+ W^ 










226 










1.7 


HSE" 










221 










3.4 


SIC-LSD^ 










214 










1.2 



M 



JPL 



** Reference [41i] . Reference |l(|, '^ Reference 



19 1, '^ Reference 1.18] j*^ Reference ij], -^Reference 



15|. 



derivative of the bulk modulus (B), all our calculated values are approximately equal to 
4.5. Overall, comparing with the experimental data and recent theoretical results, the ac- 
curacy of our atomic-structure prediction for AFM Pu02 is quite satisfactory by tuning the 
effective Hubbard parameter f/ in a range of 3-4 eV within the LDA/GGA+ U approaches, 
which supplies the safeguard for our following study of electronic structure and mechanical 
properties of Pu02. 

The total density of states (DOS) as well as the projected DOS for the Pu 5/ and O 
2p orbitals within LDA, GGA, LDA+ U, and GGA+ U formalisms are shown in Fig. |2l 



Clearly, our results reproduce all the features included in our previous work |16|. In partic- 
ular, we recover the main conclusion that although the pure LDA and GGA fail to depict 
the electronic structure, especially the insulating nature and the occupied-state character of 
Pu02, by tuning the effective Hubbard parameter in a reasonable range, the LDA/GGA+ U 
approaches can prominently improve upon the pure LDA/GGA calculations and, thus, can 
provide a satisfactory qualitative electronic structure description comparable with the pho- 




FIG. 2: (Color online) The total DOS for the Pu02 AFM phase computed in the (a) LDA, (b) 
LDA+C/ {U=A), (c) GGA, and (d) GGA+U (f/=4) formalisms. The projected DOSs for the Pu 
5/ and O 2p orbitals are also shown. The Fermi energy level is set at zero. 

toemission experiments [g, |7|. In this study, the insulating energy band gap (Eg) is of 1.5 eV 
at 11=4 eV within LDA/GGA+t/ approaches (see Table H]). In Table HTj the experimental 
10[ value and previous theoretical calculations |lj, |l5|, [l8|, [l9[ results are also listed for 
comparison. Obviously, the HSE calculations result in a larger Eg by ~1.9 eV 1^ and our 
calculated results are well consistent with the measured value and other theoretical results. 
The calculated amplitude of local spin moment is ~4.1 fiB (per Pu atom) for AFM Pu02 
within the two DFT+ U schemes and this value is comparable to previous LDA/PBE+ U 
(l9| results of ~3.9 fiB- 

To understand the chemical bonding characters of fiuorite Pu02, we present in Fig. [3] 
the crystal structure of its cubic unit cell and the charge density map of the (110) plane 
calculated in LDA+ U formalism with 11=4 eV for AFM phase. Evidently, the charge density 
around Pu and O ions are all near spherical distribution with slightly deformed toward the 
direction to their nearest neighboring atoms and there are clear covalent bridges between Pu 
and O ions. In order to describe the ionic/covalent character quantitatively and more clearly, 
we plot the line charge density distribution along the nearest Pu— O bonds (not shown). A 
minimum value of charge density (0.53 e/A^) along the Pu— O bonds is obtained and is 
listed in Table IIIII For comparison, we have also calculated some corresponding properties 
of AFM UO2 within the LDA+ U formalism. Parameters of the Hubbard term are taken as 



10 



(i-0.8 




FIG. 3: (Color online) Cubic unit cell for Pu02 in space group Fm^m. The larger blue spheres 
represent plutonium atoms and the smaller white O atoms. A slice of the isosurfaces of the electron 
density for AFM phase in (110) plane calculated in the LDA+f/ formalisms with t/=4 eV is also 
presented (in unit of e/A^ ) 



TABLE III: Calculated charge and volumes according to Bader partitioning as well as the A—0 
distances and correlated minimum values of charge density along the A—0 bonds for actinide 
dioxides. 



Compound 



QsiA) 



ClBiO) 
(e) 



Vb(^) 

(A3) 



Vb(0) A-0 

(A3) (A) 



Charge density n 

(e/A3) 



PuOa 


13.60 


7.20 


15.36 


11.58 


2.32 


NpOs 


12.52 


7.24 


15.69 


11.83 


2.34'^ 


UO2 


11.44 


7.28 


16.35 


12.05 


2.36 


ThOa 


9.34 


7.33 


17.67 


13.35 


2.43'' 



0.53 
0.5^ 
0.52 

0.45^ 



Reference [32], '' Reference 25J]. 



[/=4.51 eV and J=0.51 eV, which had been checked carefully by Dudarev et al |11|-|13|. In 
following study, results of UO2 are all calculated using these Hubbard parameters either for 
ground-state FmSm phase or high pressure Pnma phase. The lattice parameter ao=5.449 
A and bulk modulus i?=220.0 GPa for FmSm UO2 obtained by EOS fitting are in perfect 
agreement with results of recent LDA+U calculation 18| (ao=5.448 A and 5=218 GPa) 



11 



and experiments 4l|, |52| (ao=5.47 A and B=207 GPa). Charge analysis results of Fm3m 
UO2 are listed in Table lllll Clearly, the minimum values of charge density for PUO2 and 
JO2, comparable to that along the Np— O bonds included in our previous study of N pO ^ 



32|, are prominently larger than that along the Th— O bonds (0.45 e/A^) in Th02 [25 1. 
This indicates that the Pu— O, U— O, and Np— O bonds have stronger covalency than the 
Th— O bonds. And this conclusion is in good accordance with previous HSE study of the 
covalency in Ref. [2], in which significant orbital mixing and covalency in the intermediate 
region (Pu02-Cm02) of actinide dioxides was for the first time showed by the increasing 



5f-2p orbital energy degeneracy. Besides, the ^4-0 



yl=Pu, Np, U or Th) bond distances 



calculated in present study and our previous works 25|,l32] are listed in Table UTTl Obviously, 
the Pu— O, Np— O, and U— O bond distances are all smaller than the Th— O bond distance. 
This fact illustrates that the covalency of ^402 has tight relation with their bond distances, 
thus influences their macroscopical properties, such as hardness and elasticity. To see the 



ionicity of AO2, results from the Bader analysis 22|, |23|] are presented in Table IIIII The 
charge (Qs) enclosed within the Bader volume (Vb) is a good approximation to the total 
electronic charge of an atom. Note that although we have included the core charge in charge 
density calculations, since we do not expect variations as far as the trends are concerned, 
only the valence charge are listed. Here, same with our previous study, we apply [/=4.6 eV 
and J=0.6 eV for the Np 5/ orbitals in calculations of Np02. Based on the data in Table UTTl 
and considering the valence electron numbers of Pu, Np, U, Th, and O atoms (16, 15, 14, 
12, and 6, respectively), in our study of ^402, we find that about 2.40, 2.48, 2.56, and 2.66 
electrons transfer from each Pu, Np, U or Th atom to O atom, respectively. This indicates 
that the ionicity of AO2 shows decreasing trend with increasing Z. 

B. Mechanical properties 

1. Elastic properties 

Our calculated elastic constants, various moduli, and Poisson's ratio v of the FmSm 
PUO2 at GPa are collected in Table [TTl and those of the Pnma phase are listed in Table IIV[ 
Obviously, the Fm3m phase of PUO2 is mechanically stable due to the fact that its elastic 
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TABLE IV: Calculated elastic constants, elastic moduli, pressure derivative of the bulk modulus 

B , and Poisson's ratio v for cotunnite-type Pu02 and UO2 at GPa. Except B and v, all other 

values are in units of GPa. 

Compound Cn C22 C33 C44 C55 C^ Cu C23 C13 B B' Ba Bt B, G E v 

PUO2 355.6 327.4 336.3 36.0 81.5 96.7 140.5 199.0 141.2 219.9 5.7 598.5 669.5 720.3 73.0 197.2 0.351 

UO2 338.2 335.1 325.9 24.4 75.6 92.0 119.9 148.4 142.1 202.2 5.5 586.4 593.8 641.9 68.0 183.4 0.349 



constants satisfy the following mechanical stability criteria 36| of cubic structure: 



Cn > 0, C44 > 0, Cu > \Ci2\, (Cn + 2C^2) > 0. (8) 

As for the high-pressure Pnma crystalline phase, we have optimized the structural param- 
eters of its AFM phase at different pressures within LDA+ U formalism with the typical 
value of [/=4 eV. To avoid the Pulay stress problem, we perform the structure relaxation 
calculations at fixed volumes rather than constant pressures. Note that different from the 
structure optimization of the ground-state Fm3m phase, the coordinates of atoms and the 
cell shape of the Pnma phase are necessary to be optimized due to its internal degrees of 
freedom. After fitting the energy-volume data to the EOS, we obtain the optimized struc- 
tural lattice parameters a, b, and c for the Pnma PUO2 AFM phase at GPa to be 5.889, 
3.562, and 6.821 A, respectively, giving F=143.1 A^. This volume is smaller than the equi- 
librium volume of 154.2 A^ for the Fm3m phase. This feature implies that the Pnma phase 
will become stable under compression. The elastic constants listed in Table HVl indicate that 
the Pnma PUO2 is also mechanically stable. Actually, they satisfy the following mechanical 



stability criteria 36| for the orthorhombic structure: 



Cn > 0,C22 > 0,^33 > 0,^44 > 0,^55 > 0, Cge > 0, 

[Cn + C22 + C33 + 2(Ci2 + Ci3 + C23)] > 0, 

(Cn + C22 - 2C12) > 0, (Cn + C33 - 2C13) > 0, 

(C22 + C33 - 2C23) > 0. (9) 

One can see from Table HVl that the calculated C12, C23, and C13 are largely smaller than Cn, 
C22, and C33. Therefore, the mechanical stability criteria is easily satisfied. For comparison, 
we have also calculated the elastic properties of UO2 in its ground-state fiuorite phase and 
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high-pressure cotunnite phase within LDA+f/ formahsm. For FmSm phase, Cii=389.3 
GPa, Ci2=138.9 GPa, C44=71.3 GPa, 5=222.4 GPa, (7=89.5 GPa, E=236.8 GPa, and 
Poisson's ratio t;=0.323. Results for Pnma UO2 are presented in Table IIVI Clearly, both 
fluorite and cotunnite phases of UO2 are mechanically stable. Moreover, comparing results 
of bulk modulus B, shear modulus C, Young's modulus E, and Poisson's ratio v for fluorite 
phase and cotunnite phase, they are almost equal to each other for PUO2. Particularly, 
the bulk modulus B is only smaller by about 5 GPa for Pnma phase compared to that of 
Fm3m phase. As for UO2, bulk modulus of cotunnite phase is smaller by approximately 9%, 
shear modulus and Young's modulus about 23%, compared to those of fluorite structure. 



Besides, in our previous study of Th02 |25| we find that the bulk modulus, shear modulus, 
and Young's modulus of cotunnite Th02 are all smaller by approximately 25% compared to 
those of fluorite Th02. Therefore, after comparing the bulk modulus B for the two phases of 
Th02, UO2, and PUO2, we find that the difference is decreasing along with increasing of 5/ 
electrons. This trend is understandable. Although the softening in bulk modulus upon phase 
transition for these three systems is similar, the increase of 5/ electrons of actinide metal 
atoms will lead to more covalency and thus reduce the bulk modulus difference between the 
two phases of the actinide dioxides across the series. To study the anisotropy of the linear 
bulk modulus for Pnma PUO2 and UO2, we calculate the directional bulk modulus along 



the a, b and c axis by the following equations 40 1 



-^« = T"^ ri5^^b = -^,Bc = -^, (10) 

1 + a + p a p 

where A=Cii+2Ci2a+C22a^+2C*i3/3+C*33/3^+2C23a;/3. For orthorhombic crystals, 

_ [Cn — Ci2)(C33 — C13) — (C23 — Ci3)(Cii — C13) , , 

1^*33 — Ci3)(C22 — C12) — (Ci3 — C23)(Ci2 — C23) 

and 

n _ (C22 — Ci2){Cii — C13) — (Cji — Ci2)(C23 — C12) , s 

{C22 — CujyCss — C12) — (C12 — C23)(C'i3 — C23) 
are defined as the relative change of the b and c axis as a function of the deformation of 
the a axis. After calculation, we obtain values of a and /3 to be 0.8939 (0.9875) and 0.8309 
(0.9135) for Pu02 (UO2), respectively. Results of Ba, Bf,, and Be are presented in Table 
IIV[ Clearly, results of the directional bulk moduli show that both Pnma PUO2 and UO2 are 
easily compressed along a axis at GPa. The longest axis c is the hardest axis for these two 
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TABLE V: Calculated density (in g/cm'^), transverse (vt), longitudinal {vi) and average (vm) 
sound velocities (in m/s) derived from bulk and shear modulus, and Debye temperature (in K) for 



actinide dioxides. Resu^ 
our previous studies 



251, 



ts of Np02 and Th02 are calculated by using the elastic data included in 



32| 



Compound 



Phase 



Vt 



vi 



Od 



PuOa 


Fm3m 


11.892 


2477.6 


5206.3 


2787.0 


354.5 


PuOa 


Pnma 


12.812 


2387.0 


4976.0 


2683.9 


350.0 


NpOs 


Fm3m 


11.347 


2835.0 


5566.5 


3176.8 


401.2 


UO2 


Fm3>m 


11.084 


2841.8 


5552.7 


3183.4 


398.1 


UO2 


Pnma 


11.957 


2384.2 


4948.5 


2680.2 


343.7 


ThOa 


Fm3>m 


9.880 


2969.1 


5575.5 


3317.3 


402.6 


ThOa 


Pnma 


10.505 


2504.6 


4738.4 


2799.8 


346.8 



actinide dioxides. Directional bulk moduli of PuOa are all bigger than the corresponding 
values of UO2. Moreover, using results of elastic constants included in previous study |25| . 
the directional bulk moduli Ba-, Bf,, and B^. of Pnma ThOa are calculated to be 528.3, 406.3, 
and 415.8 GPa, respectively. This illustrates that, in contrary to Pnma PuOa and UO2, the 
Pnma ThOa is relatively harder to be compressed along a axis compared to other two axes. 
And all three directional bulk moduli of Pnma Th02 are apparently smaller than that of 
Pnma PUO2 and UO2, which indicates relative weaker covalency of ThOa compared with 
PuOa and UO2 in their high pressure phase. 



2. Debye temperature 

In study of the sound velocities and Debye temperature, our calculated results of AFM 
PuOa and AFM UO2 in their fluorite phase and high pressure cotunnite phase are presented 
in Table |Vl For comparison, results of Np02 and Th02 calculated by using the elastic data 



included in our previous studies 25|,l32] are also listed. As seen from Table IVl in their FmSm 



structure, Debye temperature and sound velocity of UO2, Np02, and ThOa are higher than 
that of PuOa. This is interestingly associated with the fact that Debye temperature (Od) 
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and Vickers hardness (H) of materials follow an empirical relationship 42| : 



(13) 



where M is molar mass and Q is molecular volume. Although Fm3m PUO2 has close value 
of hardness compared with UO2 and Np02, it has relatively smaller molecular volume and 
larger molar mass, as a consequence, has a lower Debye temperature than that of UO2 and 
Np02. As for the Prima structure, Debye temperature and sound velocity of PUO2, UO2, 
and Th02 are lower than those of their FmSm structure. This is since that the Pnma PUO2, 
UO2, and Th02 have smaller values of volume, bulk and shear moduli compared to their 
ground-state FmSrn structure. 



3. Hardness 

Hardness is also one fundamental physical quantity when considering the phase stabil- 
ity and mechanical properties. According to the hardness conditions 43|, the hardness is 
closely related to interatomic distances, number of nearest neighbors, directional bonding, 
anisotropy, and the indenter orientation in the structures. To date there is still no available 
calculation method involving hardness anisotropy in different dimensions in the literature. In 
spite of that, however, recently a semiempirical approach of hardness was raised by Simunek 



and Vackar 



44 



45| in terms of the atomistic bond strength properties. This approach has 



been successfully tested by the authors on the more than 30 binary structures with zinc 
blende, cubic fiuorite, rock salt crystals, and also for highly covalent crystals 4^, |45| . There 
is no need for all those high-symmetry structures to consider the anisotropy. Also, this 
method has been applied to the crystals involving covalent and ionic bonding characters, 
such as in the compounds of transition metal and light atoms, and generalized to the com- 
plex structure with more than two different bond strengths 4J, |45l]. Moreover, in study of 
covalent crystals, the results of the method raised by Simunek and Vackar 4J, |45| agree well 
with those of another hardness method of Gao and co-workers J46| . Therefore, the hardness 
of optimized cubic fiuorite structures of actinide dioxides can be calculated by the method 



of Simunek and Vackar 



44 



451]. In the case of two atoms 1 and 2 forming one bond of 



strength S12 in a unit cell of volume Q, the expression for hardness has the form 4^ 

H={C/Q)h2Sue~''^\ (14) 
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where 

si2 = \/(eie2)/(rain2(ii2),ei = Zi/vi (15) 

and 

h = C-^r = 1 - [A/(^I^/(ei + e,)r (16) 

ei + 62 

are the strength and ionicity of the chemical bond, respectively, and di2 is the interatomic 
distance; C=1550 and a=A are constants. The radius Tj is chosen to make sure that the 
sphere centered at atoms i in a crystal contains exactly the valence electronic charge Zi. 
For fluorite structure Pu02, &i2=32 counts the interatomic bonds between atoms Pu (1) 
and O (2) in the unit cell, ni=8 and n2=4 are coordination numbers of atom Pu and O, 
respectively, ri=1.698 (A) and r2=1.005 (A) are the atomic radii for Pu and O atoms, 
respectively, Zi = 16 and ^2=6 are valence charge for Pu and O atoms, respectively, c?i2=2.32 
(A) is the interatomic distance, and fi=154.16 (A'^) is the volume of unit cell. As for 
fluorite UO2, &i2=32, ni=8, n2=4, ri = 1.737 (A), r2=1.003 (A), Zi=U, ^2=6, c?i2=2.36 (A), 
^=161.82 (A3). Using Eqs. (14)-(16), we obtain Si2=0.1010 and /2=0.0503 for PUO2 and 
Si2=0.0919 and /2=0.0219 for UO2. Therefore, the hardness of PUO2 at its ground-state 
fluorite structure can be described as H=2Q.Q (GPa) and for UO2 H=25.8 (GPa). Clearly, 
the hardness of PUO2, almost equal to the hardness of Np02 (26.5 GPa) [32], is slightly 



larger than that o 
ThOs (22.4 GPa) 



UO2. Besides, these three values of hardness are all larger than that of 



25|. 



4- Theoretical tensile strength 

Although many efforts have been paid on PUO2, little is known on its theoretical tensile 
strength. The ideal strength of materials is the stress that is required to force deformation 
or fracture at the elastic instability. Although the strength of a real crystal can be changed 
by the existing cracks, dislocations, grain boundaries, and other microstructural features, its 
theoretical value can never be raised, i.e., the theoretical strength sets an upper bound on the 



attainable stress. Here, we employ a first-principles computational tensile test (FPCTT) 47 1 
to calculate the stress-strain relationship and obtain the ideal tensile strength by deforming 
the Pu02 crystals to failure. The anisotropy of the tensile strength is tested by pulling the 
initial fluorite structure along the [001], [110], and [111] directions. As shown in Fig. HI 
three geometric structures are constructed to investigate the tensile strengthes in the three 
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FIG. 4: (Color online) Schematic illustration of tension along (a) [001], (b) [110], and (c) [111] 
orientations. The AFM order is indicated by white arrows attached on Pu atoms. 
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FIG. 5: (Color online) Dependence of the (a) total energy (per formula unit), (b) stress, and (c) 
spin moments on tensile strain for AFM Pu02 in the [001], [110], and [111] directions. 

typical crystallographic directions: Hl^a) shows a general fluorite structure of PUO2 with four 
Pu and eight O atoms; \Mh) a body-centered tetragonal (bet) unit cell with two Pu and 
four O; and 111(c) a orthorhombic unit cell with six Pu and twelve O. In FPCTT, the tensile 



^ iQ^total 



where e. 



:/3 IS 



stress is calculated according to the Nielsen-Martin scheme [48 1 
the strain tensor (a, /3=1,2,3) and Q is the volume at the given tensile strain. Tensile process 
along the [001], [110], and [HI] directions is implemented by increasing the lattice constants 
along these three orientations step by step. At each step, the structure is fully relaxed until 
all other five stress components vanished except that in the tensile direction. 

The calculated total energy, stress, and spin moments as functions of uniaxial tensile 
strain for AFM PUO2 in the [001], [110], and [HI] directions are shown in Fig. El Clearly, 
all three energy-strain curves increase with increasing tensile strain, but one can easily find 
the inflexions by performing differentiations. Actually, at strains of 0.36, 0.22, and 0.18, 
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TABLE VI: Calculated stress maxima and the corresponding strains in the tensile process. 



Direction 



Stress (GPa) 



Strain 



[001] 
[110] 

[111] 



81.2 
28.3 
16.8 



0.36 
0.22 
0.18 
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FIG. 6: (Color online) Dependence of the lattice parameters on tensile strain for AFM Pu02 in 

the (a) [001], (b) [110], and (c) [111] directions. 

the stresses reach maxima 81.2, 28.3, and 16.8 GPa for pulhng in the [001], [110], and [111] 
directions, respectively. For clear comparison, all these maximum stresses (i.e., the theo- 
retical tensile strengthes in the three typical crystalline orientations) and the corresponding 
strains are listed in Table I VII Our results indicate that the [001] direction is the strongest 
tensile direction and [HI] the weakest. In fact, there are eight Pu— O bonds per formula 
unit for fluorite Pu02. The angle of all eight bonds with respect to the pulling direction 
is 45° in [001] direction. However, in [110] direction only four bonds make an angle of 45° 
with the pulling direction. Four other bonds are vertical to the pulling direction. In [111] 
direction, two bonds are parallel to the pulling direction and six others make an angle of 
about 19.5° with the pulling direction. The bonds vertical to the pulling direction have 
no contributions on the tensile strength and the bonds parallel to the pulling direction are 
easy to fracture under tensile deformation. Therefore, the fact that the tensile strength 
along the [001] direction is stronger than that along other two directions is understandable. 
Besides, we note that the stress in [110] direction experiences an abrupt decrease process 
after strain up to 0.24. This is due to the fact that the corresponding four Pu— O bonds 
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(make an angle of 45° with the puUing direction) have been pulled to fracture. The fracture 
behaviors have been clarified by plotting valence electron charge density maps (not shown). 
Under the same strain, the abrupt increase of spin moment can be clearly seen [Fig. E](c)]. 
While the spin moments in [110] and [111] directions only increase from 4.13 to 4.23 and 
4.33 fiB, respectively, the spin moments in [001] direction increase up to 5.27 ^b at the end 
of tensile deformation. In addition, the evolutions of the lattice parameters with strain in 
all three tensile processes are presented in Fig. [61 One can see from Fig. |6] that along with 
the increase of the lattice parameter in the pulling direction, other two lattice parameters 
vertical to the pulling direction are decrease smoothly. In [001] direction, the evolutions 
of the lattice parameters along [100] and [010] directions are absolutely same due to the 
structural symmetry. For all three tensile deformations, no structural transition has been 
observed in our present FPCTT study. 

C. Phonon dispersion and thermodynamic properties of fluorite Pu02 

To our knowledge, no experimental phonon frequency results have been published for 
Pu02. In 2008, Yin and Savrasov [21] successfully obtained the phonon dispersions of both 
UO2 and PUO2 by employing the LDA+DMFT scheme. They found that the dispersive 
longitudinal optical (LO) modes do not participate much in the heat transfer due to their 
large anharmonicity and only longitudinal acoustic (LA) modes having large phonon group 



49| investigated the thermo- 



velocities are efficient heat carriers. In 2009, Minamoto et al. 
dynamic properties of PUO2 based on their calculated phonon dispersion within the pure 
GGA scheme. In present work, we have calculated the Born effective charges Z* and di- 
electric constants £00 of PUO2 before phonon dispersion calculation due to their critical 
importance to correct the LO-TO splitting. For fluorite PUO2, the effective charge tensors 
for both Pu and O are isotropic because of their position symmetry. After calculation, the 
Born effective charges of Pu and O ions for AFM PUO2 are found to be Zp^=5.12 and 
Zq=— 2.56, respectively, within LDA+f/ formalism with the choice of [/=4.0 eV. In addi- 
tion, the macroscopic static dielectric tensor is also isotropic and our computed value of 
dielectric constant ^oo is 5.95 for the AFM phase. As for phonon dispersion, the Hellmann- 
Feynman theorem and the direct method ^50|] are employed to calculate the phonon curves 
along T—X—K—T—L—X—W—L directions in the BZ, together with the phonon density 
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FIG. 7: Phonon dispersion curves (left panel) and corresponding PDOS (right panel) for AFM 
Pu02 calculated within LDA+ U formalism with {7=4 eV. 

of states. We use 2x2x2 fee supereell eontaining 96 atoms and 3x3x3 Monkhorst-Paek 
A;-point mesh for the BZ integration. In order to calculate the Hellmann-Feynman forces, 
we displace four atoms (two Pu and two O atoms) from their equilibrium positions and the 
amplitude of all the displacements is 0.02 A. 

The calculated phonon dispersion curves are displayed in Fig. [7] for the AFM phase of 
Pu02. In the fiuorite primitive cell, there are only three atoms (one Pu and two O atoms). 
Therefore, nine phonon modes exist in the dispersion relations. One can see that the LO-TO 
splitting at F point becomes evident by the inclusion of polarization effects. In addition, 
due to the fact that plutonium atom is heavier than oxygen atom, the vibration frequency of 
plutonium atom is lower than that of oxygen atom. As a consequence, the phonon density 
of states for Pu02 can be viewed as two parts. One is the part lower than 6.4 THz where the 
main contribution comes from the plutonium sublattice, while the other part higher than 
6.4 THz are dominated by the dynamics of the light oxygen atoms. 

Thermodynamic properties of AFM Pu02 calculated within LDA+ U formalism with 
the choice of {7=4.0 eV are determined by phonon calculation using the quasiharmonic 
approximation 5l|], under which the Gibbs free energy G(T,P) is written as 



G(T, P) = F{T, V) + PV. 



(17) 



Here, F(T, V) is the Helmholtz free energy at temperature T and volume V and can be 
expressed as 

F{T, V) = E{V) + Fp,(T, V) + F,KT, V), (18) 
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FIG. 8: Dependence of the free energy F(T, V) on crystal lattice parameter a for a number of 
selected temperatures for AFM Pu02 calculated within LDA+ U formalism with U=A eV. 

where E(V) is the ground-state total energy, Fph (T, V) is the phonon free energy and 
Fei (T, V) is the thermal electronic contribution. In present study, we focus only on the 
contribution of atom vibrations. The Fph (T, V) can be calculated by 

FpuiT, V) = ksT j^ g{uj) In 2 sinh [—^ 

where uj=uj{V) represents the volume-dependent phonon frequencies and g{uj) is the phonon 
DOS. 

In calculation of -F(T, V), the ground-state total energy and phonon free energy should 
be calculated by constructing several 2x2x2 fee supercells. The temperature T appears 
in F{T, V) via the phonon term only. Calculated free energy F{T^ V) curves of PUO2 for 
temperature ranging from up to 1500 K are shown in Fig. |H1 The locus of the equilibrium 
lattice parameters at different temperature T are also presented. The equilibrium volume 
V(T) and the bulk modulus B(T) are obtained by EOS fitting. Figure [9] shows the tem- 
perature dependence of the relative lattice parameter and the bulk modulus. Experimental 
results J52| are also plotted. Clearly, the lattice parameter increases approximately in the 
same ratio as that in experiment [52|. The bulk modulus B(T) is predicted to decrease 
along with the increase of temperature. This kind of temperature effect is very common for 
compounds and metals. Besides, the specific heat at constant volume Cy can be directly 
calculated through 

fvjj 



a 



'0 



knT 



exp k^T 






(e^P aSt - 1)' 



du, 



(20) 



22 



0.012 
0.009 
0.006 
0.003 



s 0.000 



(a) 



^■'■' 



.n 



■ Yamashita 
A This study 



a: 






300 



600 



230 



900 1200 

Temperature (K) 



1500 




600 900 

Temperature (K) 



1500 



FIG. 9: (Color online) Temperature dependence of (a) relative lattice parameter [a{T) — 
a(300)]/a(300), where a(300) denotes the lattice parameter at 300 K, and (b) bulk modulus B(T) 
of Pu02. Experimental results from [52f] are also shown in panel (a). 
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FIG. 10: (Color online) Calculated specific heat at constant 



(Cy). For comparison, previous experimental results from 49l. 15314551] are also shown 



pressure 



{Cp) and at constant volume 



while the specific heat at constant pressure Cp can be evaluated by the thermodynamic 
relationship Cp — Cy = ay{T)B{T)V{T)T, where the isobaric thermal expansion coefficient 
can be calculated according to the formula av{T) = y (|^)p- Calculated Cp and Cv are 
presented in Fig. [TOl Clearly, the calculated Cy is in good agreement with experiment 49] 
in all investigating temperature domain, while the predicted Cp only accords well with the 
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FIG. 11: (Color online) Comparison of relative energy vs the cell volume for AFM Pu02 in Fm3m 
and Pnma phases. The total energy of Pnma phase at GPa is set as zero. A phase transition at 
24.3 GPa is predicted by the slope of the common tangent rule, as shown in the inset. 



corresponding experimental data 53j below 300 K due to the intrinsic fact that near zero 
temperature Cp and C y appro ach to the same value. The disagreement of Cp between 
theory and experiments 5J, ISSj in high temperature domain is believed to mainly originate 
from the quasiharmonic approximation we use. 



D. High pressure behavior of Pu02 

In the following, we will focus our attention on the behavior of plutonium dioxide un- 
der hydrostatic compression. Two experimentally established structures, FmSm and Pnma 
phases, are investigated in detail. The relative energies (per unit cell) of the two phases 
at different volumes are calculated and shown in Fig. [TTl Obviously, the Fm3m phase is 
stable under ambient pressure while under high pressure the Pnma phase becomes stable. 
According to the rule of common tangent of two energy curves, a phase transition at 24.3 
GPa is predicted by the slope shown in the inset of Fig. [HI Besides, we also determine the 
phase transition pressure by comparing the Gibbs free energy as a function of pressure. At 
K, the Gibbs free energy is equal to enthalpy H, expressed as H=E+PV. The crossing 
between the two enthalpy curves (not shown) also gives phase transition pressure of 24.3 
GPa, which is fully consistent with above result in terms of the common tangent rule. This 
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TABLE VII: Calculated structural parameters (in A), pressure (in GPa), and relative energy (in 
eV) with respect to the total energy at GPa of the Pnma Pu02 at two constant volumes from 



experiment 



3]. For comparison, the experimental values are also listed. 



Magnetism 


a 


b 


c 


Pressure 


Relative Energy 


NM 


6.502 


3.165 


6.087 


82.9 


32.87 


FM 


7.478 


3.088 


5.424 


124.2 


31.80 


AFM 


5.585 


3.410 


6.577 


38.6 


1.90 


Expt. 


5.64 


3.38 


6.57 


39 




Expt. 


5.62 


3.44 


6.49 


49 




AFM 


5.492 


3.398 


6.539 


48.8 


2.80 



value is smaller by ~15 GPa when compared to the experiment measurement by Dancausse 
et al 8| . We notice that Li et al 56J predicted a transition pressure of about 45 GPa em- 
ploying the full-potential linear-muffin-tin-orbital (FPLMTO) method. However, they only 
considered the NM and FM phases in their calculations of the Fm3m and Pnma Pu02. No 
results were calculated for the AFM phase. In our present study, we first optimize the Pnma 
Pu02 at constant volume (about 125.2 A'^) of the experimental value at 39 GPa |8|. Our 
calculated structural parameters of the NM, FM, and AFM phases are tabulated in Table 
IVlIi Obviously, the AFM phase is the most stable phase and its structural parameters are 
well consistent with experiment. As seen in Table IVIIl results of AFM phase at 49 GPa is 
also consistent with experiment. Therefore, we will only consider the results of AFM phase 
in the following discussion. In addition, obvious discrepancy between two works |8|, |4l| from 



the same experimental group involving the transition pressure of Th02 was noticed in our 
previous study [25]. While the transition was reported firstly to start at 40 GPa 8|], later, 
through improving experimental measurement techni que, which can capture the two phases 
cohabitation zone during transition process, Idiri et al 4]| observed that the transition really 
begins around 33 GPa. Our previous study ^^ predicted that the phase transition of Th02 
started at around 26.5 GPa. Idiri et al 411 also stated that the bulk modulus of PuOo were 

n 

argely overestimated in their previous work [8]. Thus, we believe that their former report 
Sj also overestimated the transition pressure of Pu02 in some deg ree. So the predicted 
transition pressure of Pu02, which is very close to that of Th02 [251], is understandable. 
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FIG. 12: (Color online) Partial density of states (PDOS) for (a) Fm2>m phase and (b) Pnma phase 
at around 25 GPa. The Fermi energy level is zero. The energy gaps of Fm3m phase and Pnma 
phase are 1.51 eV and 1.65 eV, respectively. 

Figure [I2] compares the PDOS of the Fm3m and Pnma phases of Pu02 at a pressure of 
around 25 GPa, close to the transition pressure. One can see evident increase of the band gap 
from 1.51 eV in the fluorite phase to 1.65 eV in the cotunnite phase. In study of UO2, Geng 
et al 2J] also predicted this similar behavior. However, our previous study of Th02 25 1 
found that the band gap is reduced from FmSm phase to Pnma phase. The reason is simply 
that Th02 is a charge-transfer insulator, while PUO2 and UO2 are the Mott-type insulators. 
From Fig [121 one can see that while O 2s and Pu 6p states expand in the low bands, O 2p and 
Pu 5f /6d states are mainly featured near the Fermi level and have prominent hybridization 
characters in a width of 6.3 (5.9) eV for Fm3m [Pnma) phase. There is no evident difference 
between the two phases in their 5/ electronic localization behavior. Therefore, our study of 
Pu02 supports our previous 25| viewpoint: the phenomenon of volume collapse during high- 
pressure phase transition of the actinide dioxides is mainly attributed to the ionic (instead 
of electronic) response to the external compression. 

The relative volume V/Vq evolution with pressure for PUO2 in both Fm3m and Pnma 
phases are shown in Fig [131 For comparison, the experimental [8l] data are also shown in the 
figure. Clearly, our calculated P-V equation of state is well consistent with the experimental 
measurement for the two phases of PUO2. Specially, at the calculated transition pressure 
(24.3 GPa), our result in Fig. [T3] gives that the volume collapse upon phase transition is 
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FIG. 13: (Color online) Calculated compression curves of Pu02 compared with experimental mea- 
surements. The volume collapses at our predicted phase transition point 24.3 GPa and experimental 
81] phase transition pressure 40 GPa are labeled. 
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FIG. 14: (Color online) Pressure behavior of the relative lattice parameters of the Pnma phase, 
where the drastic change in the relative lattice constants (region between dashed lines) indicates 
an isostructural transition. Besides, the pressure behavior of the insulating band gap is also shown. 

6.8%. This value is evidently underestimated compared with the experimental data of 12% 
8|. The discrepancy between experiment and our calculation needs more experimental and 
theoretical works to examine. After phase transition, we also find an isostructural transition 
occurring between 75 and 133 GPa for the Pnma phase. This isostructural transition of 
actinide dioxides was first found in DFT calculations of UO2 by Geng et al [2j] and then 
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FIG. 15: (Color online) The total DOS for the cotunnite Pu02 AFM phase calculated at selected 
pressures within LDA+ U formalism with (a) U=0 eV, (b) [/=4 eV, and (c) U=6 eV. The projected 
DOSs for the Pu 5/ and O 2p orbitals are also shown. The Fermi energy level is set at zero. 



observed in study of Th02 25|. The pressure dependence of the three lattice parameters 
(with respect to their values at calculated transition pressure 24.3 GPa) for the Pnma phase 



m 



25|, in 



of Pu02 are plotted in Fig. [Ml Similar to the studies of UO2 and Th02 
pressure region before 75 GPa, the responses of the three relative lattice parameters to the 
compression are anisotropic in the sense that the compression of the middle axis a is most 
rapid compared to those of the long axis c and small axis b, which vary upon compression 
almost in the same tendency. When the pressure becomes higher to be between 75 and 
133 GPa, remarkably, it reveals in Fig. [H] that all the three relative lattice parameters 
undergo dramatic variations by the fact that the small axis b has a strong rebound and the 
middle a is collapsed. When the pressure is beyond 133 GPa, then the variations of the 
three relative lattice parameters become smooth and approach isotropic compression. This 
signifies a typical isostructural transition for the Pnma phase of PUO2. 

Moreover, we also present in Fig. [13] the evolution of the insulating band gap with 
pressure for Pnma phase of PUO2. Apparently, the band gap behaves smooth in pressure 
region of 25-92 GPa, then turns to decrease suddenly from about 1.75 eV to zero under 
compression from 92 to 133 GPa. This clearly indicates that the Pnma phase will occur 
a metallic transition after external pressure exceeds 133 GPa. As Fig. [TST b) shows, with 
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increasing pressure in the crossover range between 92 and 133 GPa, the 5/ electrons in the 
cotunnite phase of Pu02 are more delocahzed and the 5/ bands are largely broadened. As 
a result, the Mott-type band gap is narrowed and even blurred, which is characterized by 
finite occupancies of O 2p and Pu 5/ orbitals at the Fermi level, by the increasing kinetic 
energy of 5/ electrons. In order to see the pressure behavior of electronic structure with 
different values of Hubbard parameter U for cotunnite Pu02, we have plotted in Fig. 15 
the total DOS and PDOS for the cotunnite Pu02 AFM phase calculated at 92, 133, and 248 
GPa within LDA+f/ formalism with U=0, 4, and 6 eV. It clearly shows that (i) the pure 
LDA always produces incorrect electronic structure in the full pressure region we considered, 
and (ii) it is not imperative to adjust U as varying pressure by the fact revealed in Figs. 
15(b) and 15(c) that the electronic properties and insulator- metal transition behavior at 
high pressures show the similar character for the two choices of Hubbard parameter 11=4 
eV and f/=6 eV. Based on this observation, in the above high-pressure calculations we have 
fixed the value of f/ to be 4 eV. It should be stressed that the metallic transition is not 
unique for Pu02. Similar phenomenon has also been observed in other actinide dioxides 



24| . Besides, the variation of the local magnetic moment of plutonium atoms is almost 
same for Fm3m and Prima phases, implying that the magnetic property is insensitive to 
the structure transition in Pu02. Actually, the calculated amplitude of local spin moment 
varies from ~4.1 to ~3.8 fiB (per Pu atom) for both fiuorite and cotunnite phases in pressure 
range from to 255 GPa. No paramagnetic transition for this material has been observed 
in present study. 

IV. CONCLUSIONS 

In conclusion, the ground state properties as well as the high pressure behavior of Pu02 
were studied by means of the first-principles DFT+ U method. By choosing the Hubbard 
U parameter around 4 eV within the LDA+ U approach, the electronic structure, lattice 
parameters, and bulk modulus were calculated for both the ambient FmSm and the high- 
pressure Prima phases of Pu02 and were shown to accord well with experiments. Results for 
UO2 were also presented for comparison. Based on these results, the Pu— O and U— O bonds 
were interpreted as displaying a mixed ionic/covalent character by electronic structure anal- 
ysis. After comparing with our previous calculations of Np02 and Th02, we demonstrated 
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that the Pu— O, U— O, and Np— O bonds have stronger covalency than the Th— O bond. 
The ionicity of Th— O bond was found to be the largest among these four kinds of actinide 
dioxides. In addition, the stabihty of the two phases at zero pressure was predicted through 
calculating elastic constants and phonon dispersion. The hardness, Debye temperature, ideal 
tensile strength, and thermodynamic properties were calculated and analyzed to support the 
practical application of Pu02. We showed that the hardness of FrnSm phase is ~27 GPa 
and the Debye temperatures of Fm3m and Pnma phases are 354.5 and 350.0 K, respectively. 
For FmSm PUO2, the ideal tensile strengths are calculated within FPCTT to be 81.2, 28.3, 
and 16.8 GPa in tensile deformations along the [100], [110], and [111] directions, respec- 
tively. The volume thermal expansion and specific heat at constant volume curves are well 
consistent with available experiments. However, the discrepancy between measured and our 
calculated specific heat at constant pressure in high temperature domain is evident. This 
needs further theoretical work by including the anharmonic ionic contribution to decrease 
this kind of discrepancy. 

In studying the pressure behavior of PUO2, we showed that the FmZm^Pnma transition 
occurs at 24.3 GPa. Although this value is large smaller than the experimental report, we 
believe that our calculated result is reasonable. One reason is that the lattice parameters of 
Prima PUO2 AFM phase calculated at around 39 and 49 GPa are precisely consistent with 
experiment. Another is the fact that the experiment needs improvement as having been 
indicated in study of UO2 and Th02 [41i]. Furthermore, we extended the pressure up to 
280 GPa for the two structures of PUO2. A metallic transition at around 133 GPa and an 
isostructural transition in pressure range of 75-133 GPa were predicted for the Pnma phase. 
Also, the calculated amplitude of local spin moment only varies from ~4.1 to ~3.8 fiB (per 
Pu atom) for both fluorite and cotunnite phases in pressure range from to 255 GPa. No 
paramagnetic transition for this material has been observed. 
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